Superconductivity in doped cubic silicon: an ab initio study 
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We study within a first-principle approach the band structure, vibrational modes and electron- 
phonon coupling in boron, aluminum and phosphorus doped silicon in the diamond phase. Our 
results provide evidences that the recently discovered superconducting transition in boron doped 
cubic silicon can be explained within a standard phonon-mediated mechanism. The importance of 
lattice compression and dopant related stretching modes are emphasized. We find that Tc can be 
increased by one order of magnitude by adopting aluminum doping instead of boron. 
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The experimental discovery by Bustarret and cowork- 
ers Q of a superconducting (SC) transition in heavily 
boron-doped silicon in the diamond phase (labeled c-Si in 
what follows) concluded a long history of research on such 
a transition in silicon-based systems and in doped semi- 
conductors in general. Until recently, a SC behavior 
in silicon had only been observed in high-pressure metal- 
lic phases such as the hexagonal or /3-tin structures 0] or 
low-pressure cage-like doped clathrate phases. This 
latter structure, which becomes superconducting upon 
heavy doping, bears much similarities with the diamond 
phase, as it is a semiconducting sp 3 network, but with a 
band gap in the visible range. [6j 

Besides doped sp 3 silicon clathrates, superconductiv- 
ity in c-Si was recently made likely by the discovery of 
a SC transition in highly boron-doped carbon diamond. 
7 1 Tunnelling spectroscopy, Q reflectivity measurements 
9| and first-principles studies within the density func- 
tional theory (DFT), performed in the virtual crystal ap- 
roximation (VCA) [l(J [ll|, [l2[ or a supercell approach, 
14, [l5| strongly suggested that the transition was 
phonon mediated. Further, the study of the SC origin in 
Si-clathrates and carbon diamond led theorists to predict 
the SC transition in doped c-Si within a crude rigid-band 
model for doping [B| and a more accurate VCA treatment 
with emphasis on boron doping. 
In this paper, we study by means of ab initio simula- 
tions within a supercell approach the electronic, vibra- 
tional and electron-phonon coupling properties of boron, 
aluminum p-doped and phosphorus ra-doped c-Si. We 
find that a standard phonon-mediated BCS approach can 
account for the experimental transition temperature ob- 
served at high boron content. We predict further that 
aluminum doping would allow to increase Tc by one or- 
der of magnitude, an effect ascribed in particular to the 
negative effect of lattice compression on Tc in the case 
of boron doping. 

The calculations are performed within a planewave 
pseudopotential implementation (IH ] of the DFT using 
the PBE functional [17| for exchange and correlation. Ul- 
trasoft pseudopotentials are used with a 20 Ry (160 Ry) 



cutoff for the expansion of the wavefunctions (charge den- 
sity), increased to 25 Ry (200 Ry) in the case of boron- 
doping. A (2x2x2) supercell containing 16 atoms is built 
with one Si atom replaced by an impurity, leading to a 
^6.25% doping concentration, in the range of the esti- 
mated 5.7-8.4 % experimental value [l| for the supercon- 
ducting boron-doped samples. The Brillouin zone (BZ) 
is sampled by a (5x5x5) k-point grid for structural relax- 
ation and the calculation of the phonon modes which are 
obtained on a (4x4x4) q-grid (k and q will refer to elec- 
tron and phonon momentum respectively). We adopt 
a much finer (10x10x10) k-point sampling for calculat- 
ing the q-dependent electron-phonon coupling constants 
A(q). A tetrahedron extrapolation technique is used to 
accelerate the summation over phonon modes in the cal- 
culation of the coupling constant A. 

We first minimize the energy of the system with re- 
spect to the cell size and atomic positions. As expected, 
the largest relaxation occurs for boron doping with an 
important contraction of the Si-B bonds (~ 2.1 A), lead- 
ing to a ~ 1.9% contraction of unit cell lattice parameter. 
Experimentally, as doped samples are constrained in the 
xy plane (parallel to the surface) by the undoped sub- 
strate, the relaxation was shown to be anisotropic with 
a t zz contraction along the z-axis ranging from 2.5% to 
3.7% from XRD analysis. Using a Poisson's ratio v of 
0.28, this yields an experimental averaged lattice com- 
pression e av = e zz (l — v)/(l + v) in the 1.4%-2.1% range. 
Our theoretical value clearly falls within these estimates. 
This lattice contraction will be shown to have important 
consequences on Tc- 

We plot in Fig.QJa) the band structure of the B:Si cell 
that we compare to that of undoped silicon at the same 
lattice parameter (dashed lines). The Fermi level lies 0.55 
eV below the top of the valence bands (VBM), showing 
that in this limit of large doping, Si:B is a degenerate 
semiconductor. Clearly, close to the Fermi level, the band 
structures of doped and pristine silicon arc very similar. 
In the case of p-doped aluminum samples, the Fermi level 
falls 0.45 eV below the VBM (Fig. [lb). 

We now turn to the vibrational properties and plot in 
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FIG. 1: Electronic DFT band structures of (a) boron-, and 
(b) aluminum-doped c-Si in the (2x2x2) cell. The Fermi levels 
are indicated by horizontal lines. Dashed lines are for the 
undoped c-Si. 



Fig. [2] the phonon density of states (ph-DOS). As 
compared to the 501 cm -1 theoretical frequency for the 
zone-center optical modes in c-Si in the present PBE ap- 
proach (520 cm -1 experimentally), one observe a ~ 40 
cm -1 softening in B-doped and Al-doped samples. In the 
case of boron-doping, we find that the effect of lattice 
parameter reduction shifts the zone-center optical modes 
of pristine c-Si to higher frequencies by 20 cm . As a 
result, the effect of B-doping really amounts to a ^60 
cm -1 softening. Further, the strong reduction of the B- 
Si bonds yields a well separated B-Si (3-fold) stretching 
mode located ^100 cm -1 above the optical modes contin- 
uum, a signature in excellent agreement with the Raman 
peak at 590-600 cm -1 observed experimentally in the su- 
perconducting B-doped samples. [l| This mode has been 
already observed experimentally in non-superconducting 
samples with a large frequency shift under uniaxial pres- 
sure related to the strong local stress around the dopant. 

EH 

The Eliashbcrg function a 2 F(u>), which allows to an- 
alyze which phonons contribute to the electron-phonon 
coupling constant A, is represented in Fig. [2ja-b) (grey 
thick line) together with the ph-DOS (black line). 
a 2 F(u>) can be calculated from the knowledge of the 
electron-phonon coupling matrix elements, [2fj| namely: 



ations in (|g q i/| 2 ) comes from the two (^-functions in Eq. 
2. By defining the nesting factor n(q), that can be ob- 
tained by setting <? q „ to unity in Eq. 2, we find that 
(\gq V \ 2 )/n(q) is smoothly varying (as a function of q) 
and well behaved at zone-center. From the explicit cal- 
culation of such a ratio on the (4x4x4) q-grid, we ex- 
trapolate (|<7 q „| 2 )/n(q) using tetrahedron techniques on 
a much finer (10xl0xl0)-grid in order to obtained a well 
converged a 2 F(uj) function. [2l[ 

The integration of the Eliashberg function yields the 
coupling constant A = 2 J du> a 2 F(ui)/uj with values of 
0.26 and 0.36 for the B-doped and Al-doped systems re- 
spectively. Our A=0.26 for B-doped c-Si is close 
to the 0.30 value found in Ref. [13 using the VCA for 
a slightly smaller 5% doping level. In the present super- 
cell approach, the contribution of the B-related stretching 
modes at high energy amounts to 15 % of the coupling 
constant. In the case of Al-doping, there is an enhanced 
contribution from modes located ~ 50 cm -1 below the 
softened optical modes and showing a significant weight 
onto the Al displacements. This is reminiscent of the 
possible importance of the B-rclatcd vibrational modes 
in the case of B:C samples as suggested experimentally 
and by theory. [IS HE! 
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FIG. 2: Phonon-DOS (black line) and Eliashberg functions 
(grey thick line) for (a) boron-, and (b) aluminum-doped c-Si 
in the (2x2x2) cell. 
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To accelerate the convergence with respect to phonon 
momentum q-sampling, we observe that the largest vari- 



Using the standard McMillan formula, or a modified 
relation accounting for the <5-like shape of the Eliashberg 
function (see Ref. 1 10r ) , and with an average phonon fre- 
quency Hluq set to the value that maximizes a 2 F(uj) (~ 
470 cm -1 ), one obtains for B:Si a transition temperature 
T c evolving from 0.24 K to 0.06 K for pi* in the 0.08- 
0.12 range. Such a temperature is clearly consistent with 
the set of values obtained in Ref. [H for the three super- 
conducting samples, with a largest Tc of 0.34 K for the 
sample with the highest boron content (~ 8.4 %). [l[ 
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The larger value of A in Al:Si can be partly attributed 
to a negative pressure effect. In the case of the Si 
clathrates, it was shown experimentally and by theory 
Q that Tc decreases under applied pressure with a 
dTc/de ~2 K. As a computational experiment, we study 
Al:Si at the B:Si theoretical lattice constant, that is we 
apply a 2.1% lattice contraction. We find that A de- 
creases by ~ 10%, an effect that can be ascribed in partic- 
ular to a ~6.5% density of states reduction under broad- 
ening of the bands by applied pressure. Further, the 
expression of A as a function of a 2 F(oj), and the inspec- 
tion of Eqs.(l-3), reveal a (l/w qi ,) 2 dependence of a given 
phonon mode to the coupling constant. This scaling re- 
duces by a factor ~ (590/470) 2 ~ 1.58 the contribution 
of the three optical bands which, upon introduction of 
boron, are shifted to a much higher energy, amounting 
to further reducing A by 0.03-0.04. 

Changes in A do not affect Tc directly as the aver- 
age phonon frequency prefactor u)q and the Coulomb re- 
pulsion parameter (i* may change as well. The limited 
evolution of N(Ep)/uj suggest that /i* should hardly 
change from one type of doping to another. The aver- 
age frequency loq decreases by ~40 cm -1 in the case of 
aluminum doping. However, this decrease does not over- 
come the increase in the parameter A and with fi* in the 
0.08-0.12 range, the transition temperature Tc is found 
to lie in between 2.8 K and 0.6 K, that is about one or- 
der of magnitude larger than the values obtained in the 
case of boron doping. Even though the effect of impu- 
rity randomness (beyond the supercell approach) and 
possible anharmonic corrections may change the A values 
calculated here, we expect this evolution of the coupling 
constant from boron to aluminum to be robust. 

We finally briefly turn to the case of hypothetical 
highly n-doped phosphorus samples. In the 6.5% dop- 
ing limit, the system is degenerate, but the conduction 
bands and the multi-valley structure of pristine c-Si are 
strongly modified, with in particular a large splitting of 
the Ti5 state. We find a A value of 0.30, not significantly 
different from the boron case and rather disappointing in 
view of early arguments in favor of multi- valley semicon- 
ductors. Q 

In conclusion, we have studied the structural, elec- 
tronic and elcctron-phonon coupling properties in highly- 
doped c-Si. Our results provide support for a standard 
phonon-mediated BCS-typc mechanism for the occur- 
rence of superconductivity in B-dopcd silicon. The tran- 
sition temperature is expected to increase by an order 
of magnitude upon aluminum doping instead of boron. 
It is not clear however if the present synthesis techniques 
(such as the gas immersion laser doping technique used in 
Ref. [H) may be adapted to other type of dopants and/or 
to a larger doping percentage. Conflicting phenomena, 
such as dopant segregation as in the case of boron doped 
diamond 24J may significantly complicate the search for 



larger Tc in c-Si. 
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